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Abstract 

We consider the implications of the small Fermi surface for the ab-plane mi- 
crowave absorption. The small Fermi surface results in two superconducting 
condensates. Linear combinations of these condensates correspond to the d- 
and g-wave pairings. Microwave electric field applied in the ab-plane induces 
Josephson transitions between the condensates. Dependence of the absorption 
upon direction, amplitude and frequency of the external microwave field is calcu- 
lated. Results of the calculation are compared with the observed recently novel 
non-linear phenomena in microwave absorption of YBaC^Ofs single crystals. 
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I. INTRODUCTION 



Recently a coherent Josephson response of the entire ab-plane has been observed in the 
microwave impedance of YBaC^O-js single crystals M. This observation can be considered 
as a direct indication of the two-condensate nature of the cuprate superconductivity. The 
purpose of the present paper is to describe the coherent microwave Josephson response within 
the theory of high-temperature superconductivity based on the t-J model and the assumption 
of a small Fermi surface. 

There is a controversy in literature about the shape of the Fermi surface in cuprate super- 
conductors. In the early days it was believed that it is the small Fermi surface of the doped 
Mott insulator ||. Later many of the results from photoelectron spectroscopy (PES) have been 
interpreted as favoring the large Fermi surface in agreement with the Lattinger's theorem ||]. 
On the other hand the most recent PES data [f|-[7j once more give indications of the small 
Fermi surface for underdoped samples. In the present paper we consider the doped Mott insu- 
lator scenario with the small Fermi surface consisting of hole pockets around (±7r/2, ±7r/2), see 
Fig. la. It is widely believed that the t — J model describes the main details of the doped Mott 
insulator. To fit the experimental hole dispersion one needs to extend the model introducing 
additional hopping matrix elements t', t" (see. e.g. Refs. |5|-l0fl), but basically it is the t — J 
model. Superconducting pairing induced by the spin-wave exchange in the t — J model has 
been considered in the papers 0,12]. It was demonstrated fll] that there is an infinite set of 
solutions for the superconducting gap and all the solutions have nodes along the lines (1, ±1), 
see Fig. la. (It is very convenient to use the magnetic Brillouin zone, but it can certainly be 
mapped to the full zone.) Using translation by the vector of the inverse magnetic lattice the 
picture can be reduced to two hole pockets centered around the points (±7r/2, vr/2), see Fig. la. 
The superconducting pairing is the strongest between particles from the same pocket, and the 
lowest energy solution for the superconducting gap has only one node line in each pocket. Hav- 
ing this solution in a single pocket, one can generate two solutions in the whole Brillouin zone 
taking symmetric or antisymmetric combinations between the pockets. The symmetric combi- 
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nation corresponds to the d-wave (Fig. lb), and the antisymmetric combination corresponds to 
the g-wave (Fig.lc) pairing. A possibility of generating new solutions by taking different com- 
binations between pockets was first demonstrated by Scalettar, Singh, and Zhang in the paper 
fT"3|| . The g-wave pairing in cuprate superconductors has also been considered by Buchholtz, 
Palumbo, Rainer, and Sauls [TJ|]. They discussed the most general case of two d-wave order 
parameters and two g-wave order parameters and their behavior near the surface of the super- 
conductor. The small Fermi surface case considered in the present paper is different, because 
the pairing exists only in one d-wave channel and in one g-wave channel. On the other hand, 
these two pairings are almost degenerate and therefore the new spatial scale Z 7 arises. This scale 
is much larger than the superconducting correlation length £. The energy splitting between the 
d- and g-wave solutions in the t — t' — J model has been investigated numerically in the Ref. 
[12| . The g-wave solution disappears, and only the d-wave one survives, as soon as the hole 
dispersion is degenerate along the face of the magnetic Brillouin zone. In this situation there 
are no pockets and one has a large open Fermi surface for arbitrary small hole concentration. 
However for small well separated pockets the d- and g-wave solutions are almost degenerate. In 
the pure t — J or t — t' — t" — J model the d-wave solution always has the lower energy. However 
if one extends the model including the nearest sites hole-hole Coulomb repulsion, the order can 
be inverted. The nearest sites repulsion does not influence the g-wave pairing and substantially 
suppresses the d-wave pairing. So it is possible that the real ground state has the g-wave su- 
perconducting gap. We would like to note that the g-wave ground state does not contradict the 



existing experimental data on the Josephson tunneling fL5fl . In this case the tunneling current 
pumps the g-wave into d-wave in a thin layer near the contact, giving the interference picture 



which is very close to that for a pure d-wave ground state [jig] . For the present work it does not 
matter whether the ground state has d-wave or g-wave symmetry. What is important is that 
there are two coupled superconducting condensates. Then the idea is straightforward: external 
microwave field induces Josephson transitions between these condensates. We stress that the 
condensates are not separated in coordinate space, they coexist at each point and are separated 
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only in k-space. 

The possibility of two close condensates has been suggested for conventional superconductors 



a long time ago by Suhl, Matthias and Walker [IT]]. Collective excitations in such a system 



were considered by Legget ]nj , and then Rogovin discussed the behavior of the two component 
system in the external microwave electric field [B5|. In the present work we consider the response 
of the cuprate superconductor to an external microwave electric field parallel to the ab-plane. 
Small orthorhombic distortion is not important for this process, and therefore it is neglected. 
Hence we neglect twinning and consider a single crystal. 



II. THE GINZBURG-LANDAU LAGRANGIAN AND FREE ENERGY 

According to the microscopic picture described above we should consider simultaneously the 
d- and g-wave pairings. Let us formulate an effective Ginzburg-Landau theory describing this 
situation. In the first approximation we can neglect the interaction between pockets in k-space. 
Then half of the holes belong to one pocket and the rest belong to the other pocket, and we 
should introduce two macroscopic condensates corresponding to the pockets, = l^ile*^ 1 , 
\&2 — 1^2 1 e^ 2 , where = j^l — \/iV^/2, is the number density of condensate holes. The 
effective Lagrangian of the system in an external electric field E is of the form (hereafter we 
set h = 1) 



L = E \ {%% - %%) - F 



p=l,2 



(1) 



where F is Ginzburg-Landau free energy 
1 



E 



p=l,2 \*" l p 



No/4: 



+ ^\ dV + F int, (2) 



where if is a scalar potential, N is an equilibrium number density of condensate holes, and 
-^int i s a sman interaction between pockets. We choose axes x and y to be directed at 45° with 
respect to the crystal axes a and b. For such choice the mass tensor in each pocket is diagonal 



and m" a = x,y are corresponding masses. Due to symmetry of the pockets = m| = m 



and m\ = m\ = M. In eq. (0) and hereafter we assume summation over a = x, y. Interaction 
between the pockets is described by term Fj nt in eq.(0). Following Legget |T5[ we use the 
simplest form of this interaction 

F- m t=lj m^2 + ^i%)dV, (3) 

where 7 < a is a small parameter of the interaction. For the homogeneous case 
-^int ~~ * ^""fV I I ^2 1 cos(0i — 02), where V is the total volume. The Hamiltonian corresponding 
to the Lagrangian ([!]) is the Ginzburg-Landau free energy (Q). The equilibrium values of the 
order parameters are 

,9 , ,9 , ,9 cl + I7I Nn a , . 

I 11 I 21 1 1 2fe 4 26 V ; 

The ground state phase difference = 02 — 0i is determined by the sign of 7. If 7 > then 
= 7r and the ground state pairing has g-wave symmetry. If 7 < then = and the ground 
state pairing has d-wave symmetry. 

Note that the free energy (fj) is written in the long- wavelength limit. This approxima- 
tion is valid if the wavelength of the microwave field is much larger than the superconduct- 
ing correlation length £ = (2ay/ mM^j ^ as well as the scale related to the weak coupling 
ly = {^yy/mM^j : X 3> / 7 ,£. In this limit it is convenient to use a gauge with the scalar 
potential of the external field ip = — E • r and the vector potential A = |H x r. In this gauge 
the effect of the vector potential is small and therefore it is neglected in the free energy (fj). 
Note also that there is an obvious upper limit on the applied electric field: eE^ <C A, where 
A ~ 300K is a typical value of the superconducting gap. 

III. THE EQUATIONS OF MOTION AND THE BULK JOSEPHSON EFFECT 

Lagrange equations corresponding to ([I]) are 

Aip = Sire (l^) 2 + |^ 2 | 2 - N /2) , 
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where p = 2 if p=l, and p = 1 if p=2. To be specific let us consider the case 7 < which 
corresponds to the d-wave ground state. Stiffness of the free energy (0) with respect to variation 
of |* p | is much larger than that with respect to the phases variation. Therefore we can represent 



the order parameters as * p = yiVo/4(l + 5|* p |) exp(z0 p ), where <5|* p | <C 1, but the phase P 
is arbitrary. In the linear approximation in 5|* p | the eqs.([5]) can be rewritten as 
d5\ty p \ 1 



dt 2m* V 



(2V Q p V a <5|* p | + V 2 p ) + |7| sin(0 p - <j> p ), 



dip 



V 



dt 2m p 



l — (V^|* p | - (V Q P ) 2 ) + WV 5|* P | + 2eip - | 7 | (cos(0 p - <f> p ) - 1) , (6) 



Ay? = -47reiV (5|^i| + 5|* 
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We will see below that <5|*i| + <5|* 2 | = 0, i.e. there is no charge density. Therefore the scalar 
potential is due to the external field only: <p = — E • r • sin cut = —(E x x + E y y) sino;t, where E 
and uo are the amplitude and the frequency of the microwave field. Let us introduce the phase 



2e F e 2 (El El 

which satisfies the equation 



E . r . coswt _ (2wt + sin2wt), (7) 

y u 2ui 6 \m* m p J 



— (V a 0j) 2 -2eE-r-sin(u;t). (8) 

7l p 



dt 2m 

After substitution P = p + <50 p i n t° ec L s - © we obtain equations for <50 p and <5|* p |. We 
are looking for <50 p and <5|* p | independent of coordinates and therefore the terms with spatial 
derivatives can be omitted resulting in 

±{5\* 1 \ + 6\* 2 \)=Q, 

j t (5|* 2 |- 5|*i|) = 2| 7 | sin (9) 
^(50 2 + 50 1 ) = -4| 7 |sin 2 (0/2), 

-|(<50 2 -6(h) = bN (5\* 1 \-5\* 2 \), 
at 

where = 02 — 4>i = 02 ~ 4>i + 502 — 50i. First of eqs.@ proves that there is no induced charge 
density: 5|*i| + 5|* 2 | = 0. According to eq. (0) 



0° - 0? = -A(2ut + sin2^), A = ^ (e> - Ef) (10) 

From @ and fllPf ) we find the equation for the relative phase 

(f> + cul sin (p = A(2to) 2 sin 2ujt, (11) 

where ^ = 2^J~fy\a is the threshold frequency of the phase excitation in the d-g-wave cuprates 
20|j . Note that the driving force in this equation reminds of that found in the paper [|19f| . 
Equation (|TTD does not contain absorption. In the present work we follow the simplest 



phenomenological way of just adding a "friction term" Tcj) to the left hand side of eq. (|TT|) 
rewriting it as 

+ T0 + ^sin0 = A(2uof sin2uot. (12) 

^From this equation one immediately derives that absorption is proportional to 

R oc (r0 2 ) oc Auj 3 / <j)sm(2ut)dt. (13) 
Jo 

The form of eq. ([l2|) is identical to the one describing the resistively shunted Josephson junction 
which has been used in the work ffl] for phenomenological description of experimental data, see 



also Ref. . We have derived this equation from the microscopic picture based on the small 
Fermi surface. The derivation gives us all coefficients except the "friction" which has been 
introduced phenomenologically. This approach allows to find the dependence of the absorption 
upon direction, amplitude and frequency of the external microwave field. 
Concluding this section let us calculate the electric current. 

j a = ^|^i| 2 V Q 0i + ^|^2| 2 V a 2 . (14) 

In the phase (fi p = 0° + 54> p only 0° is dependent on coordinates. Using eq. (0) one finds 

j = -M(I+ 1 Ve-cosu,*. (15) 

This is the usual time dependent Meissner current. 



IV. NUMERICAL ESTIMATIONS AND COMPARISON WITH EXPERIMENT 



Dimensionless form of eq. ( |TT| ) is 

A0 + + sin0 = a sin fir, (16) 

where r = tu 2 /T, A = u 2 /T 2 , fl = 2uT/lJq, a = 4Au 2 /uq, and the "dot" denotes the derivative 
over dimensionless time r. Equation (|16|) with A = ( "overdamped pendulum" ) has been used 
in Ref. |IJ to describe the experimental data on coherent Josephson response of Y BaCu20js 
single crystals at microwave frequency to = 27rl0 10 . The data has been fitted with Q pa 0.1, 
and we take this value as an experimental result. Unfortunately our microscopic calculation 
does not give the value of T. A reasonable fit of the experimental data [[]]] can be obtained with 
r > uj . To be specific we set T = ujq, then uj = 2uj/Q = 2tt ■ 200GHz, which is close to the 
estimation ujq ~ 2ir ■ 300GHz obtained in Ref. [[2tJ from the microscopic picture. Results of a 
numerical calculation of the normalized absorption R = {^sinfir^ /a are shown in Fig.2. All 
features of the impedance data [§]] are reproduced and the first step occurs at \a\ ~ 1.1. (The 
critical value of \a\ ~ 1 is evident from eq. fll6|) without any computations.) This allows one 
to find the critical field E c - the field at which the first step occurs. Using eq. (|10"D we find 

|a| pa 1. (17) 



2e 2 M-m - 

E z c cos (2(3) 



uju 1 mM 

Here f3 is the angle between the field and the axis x or y. Taking cos (2/9) = 1, reduced mass 
mM/ (M—m) equal to the electron mass, and the frequencies uj = 2n-10GHz, uo = 2tt-200GHz 
we find the value of the critical field: E c ~ 450mGs. The experimental critical field is H c = 
200mGs [|l[. This is in good agreement if we assume that at the surface of the sample the 
magnetic microwave field is of the order of the electric one. Equation (|P7|) determines the 
dependence of the critical field on the microwave frequency and the angle between the field and 
the x axis 



E r CX 



UJ 



cos(2/3)| v ' 

Note that | cos(2/5)| = | sin(2#)| where 9 is the angle between the electric microwave field and 

the crystal axes a or b. 
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V. CONCLUSIONS 



We have considered the scenario with the small Fermi surface consisting of hole pockets. The 
small Fermi surface together with the mechanism of the magnetic pairing result in the possibility 
of having two superconducting condensates. Energy splitting between these condensates is small 
and their linear combinations correspond to the d- and g-wave pairings. The ground state 
symmetry depends on the interplay between the magnetic pairing and the Coulomb repulsion. 
Microwave electric field applied in the ab-plane induces Josephson transitions between the 
condensates. This results in a nonlinear absorption when the entire crystal responds like a 
single Josephson junction. This could be an explanation of the effect observed recently in 
YBaCu20rs single crystals 0. The theoretical estimation for the critical microwave field 
obtained in the present work H c ~ E c ~ 450mGs is close to the experimental value. We 
predict the dependence of the critical field on the microwave frequency and orientation of the 
microwave electric field with respect to the crystal axes of the cuprate. This dependence is 
given by eq. (|l"8D . 
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FIGURE CAPTIONS 

Fig.l.a. Fermi surface in magnetic Brillouin zone which is equivalent to the two-pocket 
Fermi surface (dashed line), b. Symmetry of the d-wave pairing in momentum space, c. 
Symmetry of the g-wave pairing in momentum space. 

Fig. 2. Normalized absorption R versus dimensionless driving force a. 
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